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ABSTRACT 


\ 


Exact  solutions  for  the  reflection  and  transmission  coefficients  are 
obtained  for  two  general  electron  density  distributions:  kinked-trapesoid 
(Fig.  1)  and  exponential -homogeneous -exponential  (Fig.  5).  Normal 
incidence  into  a  stratified  plasma  slab  is  assumed  so  that  the  electromagnetic 
waves  are  propagating  parallel  to  the  free  electron  density  gradients.  A 
constant  magnetic  field  is  applied  in  the  propagation  direction.  The  solutions 
derived  on  this  basis  are  used  to  evaluate  the  effect  of  the  more  realistic 
inhomogeneous  plasma  model,  parametrically.  Specific  consideration  is 
given  to  the  analysis  of  transmission  from  a  re-entry  cone.  The  asymptotic 
expansion  of  the  exact  solution  for  the  kinked -trapezoid  profile  provides  a 
simplified  expression  for  calculation  purposes  when  L/XQ  is  larger  than  one. 
The  analysis  of  problems  in  which  L/X^  is  much  smaller  than  one  is  con¬ 
sidered  in  some  detail  in  the  final  section. 
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INTRODUCTION 

The  general  problem  of  the  propagation  of  harmonic  plane  electro* 
magnetic  waves  in  magnetoactive  plasmas  was  considered  in  Refs.  1  and  2. 

A  brief  review  of  the  fundamental  equations  and  a  discussion  of  the  problem 
of  finding  specific  analytical  solutions  were  presented.  Detailed  conclusions 
as  to  the  effect  of  the  applied  magnetic  field  on  transmission  were  based, 
primarily,  on  the  parametric  analysis  of  a  homogeneous  plasma  slab. 
Preliminary  consideration  was  also  given  to  the  more  realistic  inhomogeneous 
plasma  model. 

In  a  number  of  applications  of  current  interest,  the  plasma  contains 
free  electron  density  gradients.  For  example,  in  microwave  diagnostics  of 
ionized  gas  flows  in  shock  tubes,  plasma  tunnels,  and  fusion  machines,  the 
electromagnetic  wave  is  caused  to  propagate  from  a  dielectric  into  a  dissipa¬ 
tive  medium  through  thermal  boundary  layers  having  such  gradients.  A 
similar  situation  is  encountered  by  radar  signals  propagating  to  or  from 
re-entry  bodies.  In  this  case,  the  radar  signals  may  intercept  ionized  wakes 
also  which  are  themselves  inhomogeneous.  The  same  problem  exists  in  the 
now  classical  consideration  of  long  radio  waves  entering  the  ionosphere. 

In  each  case,  if  the  electron  density  is  slowly  varying,  a  WKB  approxi¬ 
mation  can  be  used  to  represent  the  propagation  in  the  inhomogeneous 
region.  If  the  transition  distance  between  the  dielectric  and  the  homogeneous 
region  is  small  compared  to  the  plasma  slab  thickness,  as  is  very  often  the 
case  in  the  aforementioned  laboratory  work,  then  it  is  reasonable  to  assume 
complete  homogeneity.  In  general,  however,  it  will  be  necessary  to  examine 
the  complete  inhomogeneous  plasma  problem. 


*This  term  is  used  to  conveniently  indicate  the  presence  of  an  applied 
magnetic  field. 
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The  purpose  of  this  report  is  to  obtain  the  reflection  and  transmission 
coefficients  for  several  exemplary  transition  zones  between  dielectric  and 
dissipative  gases.  A  number  of  formal  solutions  have  been  obtained^  particu¬ 
larly  in  radio  ionospheric  research.  However,  relatively  few  boundary  value 
problems  have  been  treated  in  detail.  We  will  restrict  our  attention  to  the 
case  of  normal  incidence  into  a  stratified  medium  such  that  the  electro¬ 
magnetic  waves  are  propagating  parallel  to  the  free  electron  density  gradients. 
The  magnetic  field  is  constant  and  is  applied  in  the  propagation  direction. 

The  solutions  derived  on  this  basis  are  used  to  evaluate  the  effect  of  the 
more  realistic  inhomogeneous  plasma  model,  parametrically,  both  with  and 
without  an  applied  magnetic  field.  At  the  same  time,  previously^  obtained 
homogeneous  plasma  calculations  for  typical  flight  conditions  can  be  refined 
accordingly.  It  should  be  noted  that  although  the  present  problem  is  formu¬ 
lated  in  terms  of  ionized  gas  parameters,  the  equations  and  results  are 
formally  applicable  to  a  variety  of  problems  in  other  areas  of  interest  — 
for  example  in  the  study  of  certain  electromagnetic  properties  of  solids. 
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SECTION  II 


FUNDAMENTAL  EQUATIONS 

The  electromagnetic  properties  of  a  slightly  ionised  gas  may  he 
characterised  by  means  of  &  conductivity  function  which  linearly  relates  the 
current  density  and  electric  field  strength.  This  function  can  be  determined 
from  the  equation  of  motion  for  an  average  electron  in  some  region  of  space 
in  which  is  established,  in  general,  a  steady  spatially  dependent  biasing 
magnetic  field  of  induction  : 


m|~  +  mwcv  =  e(ff  +  v  X  SQ)  ,  (1) 

where  m  is  the  electron  mass,  w  is  the  average  electron  collision  frequency, 
\T  is  the  electron  drift  velocity,  e  the  charge,  and  if  is  the  electric  field 
intensity.  The  force  exerted  by  the  magnetic  field  of  the  wave  is  neglected, 
harmonic  time  dependence  of  the  form  exp(ia>t)  for  all  time -varying  vectors 
is  assumed,  and  rationalized  mks  units  are  used. 

A  right-handed  Cartesian  coordinate  system  with  the  positive  z-axis 
vertical  is  chosen,  and  the  planes  z  *  0  and  z  =  L  constitute  abrupt  bound¬ 
aries  between  free  space  in  the  regions  z  <  0  and  z  >  L  and  the  plasma  in 
the  region  0  <  z  <  L.  The  plasma  is  assumed  to  be  stratified  so  that  its 
properties  are  functions  only  of  z.  We  consider  a  plane  wave  from  free 
space  to  be  incident  normally  at  the  abrupt  free  space-plasma  boundary  z  =  0, 
propagating  into  the  medium  with  its  wave -normal  along  the  positive  z-axis. 
Further,  we  assume  that  the  applied  magnetic  field  is  in  the  direction  of 
propagation  and  that  the  induction  is  constant.  Under  the  conditions  noted 
above  and  using  the  relation 

7  =  Nev  ,  (2) 
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where  N  is  the  number  density  of  electrons ,  we  can  solve  Eq*  (1)  lor  the  current 
density  If.  Substituting  this  expression  into  Maxwell's  equations*  we  obtain 
the  following  governing  equations  for  the  propagation  of  the  right-  and  left- 
handed  circularly  polarised  waves* 


where  n^  =  w/cq  =  2ir/XQ  is  the  free  space  wave  propagation  constant,  Xq  is 
the  free  space  wavelength,  and 


( 


and  the  vector  f  act*  tn  the  opposite  direction  to  (sine*  «  <  0) 
where  Y  *  1 7 1  ■  |e  |  Bp/me*. 

In  the  following  section*,  solution*  are  obtained  to  Eq*.  (3)  and  (4)  for 
several  boundary  value  problem*  of  interest  in  which  the  electron  density 
distribution  and,  therefore  X,  is  a  known  function  of  a .  The  collision 
frequency  is  assumed  constant  in  all  cases.  2  The  general  procedure  for 
deriving  these  solutions  and  the  specification  of  appropriate  boundary 
conditions  will  be  brought  out  in  the  separate  analyses;  note,  however,  that 
if  we  let 


i 


mr  =  (U-Y)*1  ,  m|*(U  +  Yf1  ,  (S) 

then  the  equations  for  both  the  right*-  and  the  left-hand  waves  can  be  written 
in  the  form 


F"  +  n*{  1  -  mX(i)] F  =  0 


(9) 


G  = 


(10) 


( 
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SECTION  in 


PROPAGATION  IN  INHOMOGENEOUS  PLASMAS 


A.  Kinked- Trapezoid  Electron  Density  Distribution 

In  order  to  evaluate  the  influence  of  an  inhomogeneous  plasma  region 
on  propagation,  it  is  clearly  desirable  to  be  able  to  vary  die  character  of  the 
assumed  profile  while  holding  the  remaining  physical  parameters  fixed. 
Perhaps  the  simplest  way  of  accomplishing  this  is  shown  in  Fig.  1.  *  The 
transition  regions  0  <  z  <  z^,  z^  <  z  <  L  connecting  free  space  (z  <  0,  %  >  L) 
with  a  homogeneous  plasma  (Zj  <  z  <  z^)  are  represented  by  two  straight-line 
segments.  For  simplicity  die  profile  is  assumed  to  be  symmetric,  i.  e.  , 

=  L  -  Zq  and  z^  =  L  -  z^.  By  varying  the  normalized  transition  dimension 
z^/L,  the  location  of  the  kink  Zq/zj,  and  the  value  of  X  at  the  kink  X^/Xj, 
while  holding  Xj,  Z,  Y,  and  L /\q  fixed,  we  can  isolate  the  effect  of  the  in¬ 
homogeneous  plasma  regions  on  reflection  and  transmission,  parametrically. 


Consider  the  following  boundary  value  problem  consistent  with  the 
description  of  Section  II: 


F"  +  n£[  1  -  mX(s)]  F  =  0 
mr  =  (1  -  Y  -  iZf1 

mJ  »  (1  +  Y  -  iZf1 


<n> 


*This  approach  w»i  suggested  by  Dr.  R.  Jahn  in  the  course  of  several 
discussions  of  the  general  problem. 


1 


( 


X(a)  «  0 


s  <  0 


*  XQ  s/ *0  ,  0<b<b, 


i  *  -  “0 


*  Xn  +  (Xj  -  Xn)(s  -  *n)/(»|  -  *ft) 


1 


*1<  *<  *2 


=  XQ  -  (Xj  -  Xq)(b  -  Bj)/^  -  b0) 
=  -Xq(b  -  L)/bq  , 

=  0  ,  z  >  L 


*3  <  *  <  L 


*0<  *!*i 


*2  i  *  1  *3 


(12) 


The  solutions  in  the  several  regions  are  as  follows.  For  z  <  0 


F(z)  -  c'1*  t  c<2>  .‘"O' 


(13) 


For  0  <  z  <  Zq,  Eq.  (11)  becomes 


if  mX0R\ 

+ n°\'  -  —r  - 


Introducing  the  change  of  independent  variables 


-m^)  ■ 


we  obtain  Stokes'  equation 


d2F 

dP 


=  tF 


(14) 
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and*  therefore,  a  solution  in  terms  of  Airy  functions 


FU)  »  c(3Wi(fe)  +  c(4Wi(|) 


(15) 


If  we  denote  the  index  of  refraction  at  the  point  z  *  x.  by 


K2  =  n*2(n.)  =  1  -  mX(z.)  =  1  -  mXj 


(16) 


then  in  a  similar  manner,  the  substitution 


fnOK0(*i  -  V?'!  /,  KlV“  *  *o\l 

It^iJ  l  v  '^~/J 


(17) 


results  in  the  solution,  for  <  z  <  z^t 

F(r,)  =  c(5)*i(ti)  +  c(6)#i(tl) 

For  the  remaining  regions:  Zz  —  *  £  *3*  z3  £  as  <  L»,  z  >  L, 

respectively,  the  following  solutions  are  obtained: 

F(.)  .  «<7>  ♦  «(»  >K|1  . 


F(t)  «  c(9)*i(0  +  c(10Wi(4)  , 


(18) 


t-- 


*0K0<*1  *  ‘O^T,  .(,  K?V*‘Z3\ 


(19) 


(20) 


(21) 


t 


r(v)»e,ll,4»t(»)  +  e|ll<#l(v)  ,  02) 


F(.)  .  c'l3»  a'*"0' 


Continuity  of  F(x)  and  G(z)  [or,  equivalently,  F'(z),  in  view  of 
Eq.  (10)]  at  *  =  0,  z^9  z^9  z^t  and  L  gives  rise  to  12  equations  for  the 

11  amplitudes  and,  hence,  the  two  ratios  of  interest  (for  each  wave), 

R  .  c'2>/c<‘>  and  T  =  c«13'/c»». 


( 


where  the  upper  sign  on  Q  identifies  the  sign  on  both  terms  in  the  first 
column  while  th^  lower  sign  fixes  the  sign  on  both  terms  in  the  second  row, 
and 


A+(*)  A‘(¥) 

A+(*)  A‘(lr) 

pr 

A+(»K0)  A'(»K0) 

-  P2  = 

B+(*K0)  B*(#Kq) 

B+(*)  B‘(¥) 

B+(¥)  B  *(♦) 

P,  = 

.  P,  « 

3 

A+(#K0)  A‘(*K0) 

'  4 

B+(*Kq)  B'(*K0) 

The  basic  composite  functions  in  the  solution  are,  for  any  argument  P, 


i 


A*(fJ)  =  Ai'  (P2)  ±  PAi(p2) 
B*{p)  =#i'(p2)  ±  P*i((J2)  , 


and  the  arguments  are  of  the  form, 


(27) 


(28) 


The  ratios  of  the  reflected  and  transmitted  energy  to  the  incident 
energy, 

(29) 

(30) 

( 


*R  'j(lRr|2+  lRil2) 
.T=l(lTr|2+  IT(|2) 
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**m  computed  in  a  parametric  fashion  from  the  lolvtioM  fit)  and  (id).  The 
incident  wave  (a  <  0)  ia  now  a* turned  to  be  linearly  polarised,  aay  along  the 
x-ajds  (Es  and  Hy  components  only),  so  dint  the  amplitude  is  die  same 
for  die  left*  and  right-hand  waves.  The  principal  effect  of  interest  is  that  of 
die  inhomogeneity  which  manifests  itself  through  two  independent  quantities, 
namely,  die  transition  distance  and  the  profile  details.  Hie  ratio  of  transi¬ 
tion  distance  to  slab  thickness  has  the  values  0  <  a j  /L  <  0.  5,  in  die  present 
problem,  while  16  separate  profiles  are  assumed  as  determined  by  die  loca¬ 
tion  of  the  kink  and  its  ordinate,  i.  e.  ,  Zq/z^  =  0.  2,0.  4,0.  6,0.  8,  and 
VXj  =  2.0.  4,0.  6,0.  8.  This  effect  is  demonstrated  by  plotting  the  deg¬ 

radation  in  transmitted  energy  db^.  s  10  log  jq  i^,  normalised  with  respect 
to  die  corresponding  homogeneous  plasma  slab  value  (s^/L  *  0),  versus  Zj/L 
for  the  extreme  cases  XQ/Xj  =0.2,  zQ/zj  =0.8  and  Xg/X j  •  0.  8,  Zq/zj  =  0.  2, 
as  well  as  for  the  straight  trapezoidal  profile  Xq/Xj  *  Zq/zj.  Only  numerical 
results  for  transmission  are  presented  in  this  report. 

It  is  immediately  apparent  that  the  transmission  through  an  inhomo¬ 
geneous  plasma  slab  can  be  substantially  greater  than  the  corresponding 
value  for  the  homogeneous  slab.  Both  die  transition  distance  and  the  profile 
details  can  significantly  reflect  this  effect.  Thus,  in  Fig.  2a (Z  =  1),  for 
example,  when  z^/L  =  0.  5,  transmission  is  increased  by  21  percent  over  die 
corresponding  homogeneous  plasma  result  for  Xq/X^  =  0.  8,  Zq/zj  =  0.  2  and 
by  83  percent  for  Xq/Xj  =  0.  2,  z^/zj  =  0.  8.  In  die  latter  case,  a  10  percent 
increase  is  shown  for  an  order  of  magnitude  smaller  value  of  z^/L  *  0.  05. 

A  qualitative  appreciation  for  what  is  involved  can  be  obtained  from  the  as¬ 
sumed  distribution  of  X(z)  shown  in  Fig.  1.  In  particular,  one  would  expect 
to  improve  transmission,  in  general,  when  die  plasma  is  made  less  dense 
relative  to  the  homogeneous  case.  This  is  essentially  what  is  done  when  z^/L 
is  increased  (for  fixed  values  of  Xj,  Z,  and  L/X.q)  or  when  the  assumed 
profile  is  altered  by  decreasing  Xq/X^  or  increasing  Zg/Z|.  The  accuracy  of 
a  homogeneous  plasma  slab  analysis  will  improve  as  the  distribution  of  X(s) 
approaches  the  constant  value  X  j . 
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Clearly,  transmission  also  d spends  on  the  reference  parameters  Xj, 

Z,  Y,  and  L/XQ.  For  the  same  conditions  noted  above ,  for  example ,  a 
value  of  Z  *  0.  01  results  in  increases  in  transmission  of  SI  percent  and  % 
percent,  respectively.  The  influence  of  Z  in  this  figure  is  even  more  sob* 
stantial  for  the  trapezoidal  profile.  Parametric  calculations  were  made  for 
L/\0  a  0.  5  and  1. 0  including  moderately  overdense  conditions  (i.  e.  f 
1  <  <  4),  comparable  values  of  0  <  Y  <  3  required  to  improve  transmission 

to  the  level  of  practical  interest,  and  0.  01  <  Z  <  1.  Illustrative  results  are 
shown  in  Figs.  2a  through  2k.  Estimates  can  be  made  in  any  particular 
application  with  this  series  of  graphs.  Additional  calculations  would  be  re¬ 
quired  if  the  conditions  of  interest  are  too  far  removed  from  these  values. 
This  situation  will  be  considered  further  in  the  analyses  of  Sections  III.  B  and 
IV  as  noted  in  the  following  paragraph. 

Several  curves  in  Figs.  2f  and  2g  are  incomplete  in  that  the  results 
are  not  shown  for  the  larger  values  of  z^/L,  which  calls  attention  to  the 
fact  that  the  computations  using  Eq.  (26)  met  with  numerical  difficulties  when 
the  magnitude  of  any  of  the  arguments  •,  %  #Kq,  or  PK^/Kq  became  large. 
This  difficulty  not  only  precluded  the  afore  mentioned  calculations  at 
Zj/L  =  0.  5  but  also  prevented  the  desirable  extension  of  our  parametric 
analysis  to  larger  values  of  and  Y,  as  well  as  L/Xq.  Since  highly  over- 
dense  plasmas  are  encountered  in  a  number  of  applications,  in  which  case 
values  of  Y  of  comparable  size  would  be  required  to  make  the  applied  mag¬ 
netic  field  effective,  it  is  desirable  to  consider  the  range  Y  »  1.  This 
information  is  derived  from  the  analysis  of  the  next  section. 

It  should  be  remarked  that  in  the  cases  where  at  least  one  argument 
becomes  too  large  for  our  calculation  procedure,  the  conditions  of  interest 
at  L/Xq  s  0.  5  and  1. 0  are  such  that  one  or  more  of  the  remaining  arguments 
may  still  be  relatively  small,  thereby  forestalling  any  general  asymptotic 
analysis.  An  asymptotic  expansion  of  the  exact  solution  is  derived  below, 
however,  to  provide  information  for  larger  values  of  L/Xq.  In  this  regard, 
note  that  expansions  of  any  solution  for  large  (or  small)  values  of  more  than 
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on#  independent  quantity  in  any  problem  are  meaningful  only  when  proper 
consideration  ie  given  to  their  order  of  magnitude  relative  to  each  other  a* 
well  at  relative  to  the  number  one.  This  formal  observation  becomes  of 
practical  interest  in  the  proper  use  of  the  resultant  simplified  expressions 
for  computational  purposes.  The  small  L/Kq  limit  will  be  examined  sep¬ 
arately  in  Section  IV. 

Using  the  asymptotic  representation  of  the  Airy  function,  we  obtain 


A±W  =  'XP("3p3)(a0  *  *1S’3  +  a2P 


)  ' 


B*(M  =  b*p-J  +  b*|T6  +  b*p"9  +  ■  •  , 


for  -tt  <  argP  <  tt,  where 


=  -1  ±  1 


a*  =  -(7  ±  5)/48  =  b* 


a*  =  (455  ±  385)/4608  =  -b* 

a*  =  -(95095  ±  85085)/663552  =  b* 

bj  =  1  *  1 
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The  expansion  ceases  to  be  valid  when 


0  =  ^  -  Zkir 


or  9  +  3  arg  Kc 


3ir 

:  X  *  ^ 


k  =  0,1,2, 


-1 


where  0  »  -tan’  Z/M;  and  convergence  will  be  poor  in  the  vicinity  of  these 
conditions.  The  general  asymptotic  expansion  of  (25)  and  (26)  for  arbitrary, 
large  values  of  the  parameters  becomes 


v2(SjS3  -  s351)[s453  +  +  (S2?4  -  S4S2,[S354  +  (1/64)Sj?2] 

"  8[S3S4  +  (1/64)SjS2J^  -  (y2/8)(SjS3  -  S^)2 


(31) 


T  ~ 


■y  exp(in0L) 


[S3S4  +  (l/MJSftF  -  (//64)(S1S3  -  S  jgj4 


(32) 


where 

Sj  =  V_(*)V+(®K0)tf'  -  VJ*K0)V+(*)f  + 

S2  =  W  (#K0)W+(#)f +  -  W  («)W+(*K0)«  ‘ 
s3  =  W_(*K0)V+(»)£+  +  gy  V_(*)W+(*K0>tf  ‘ 
S4  =  W_(#)V+(*K0)£  *  +  ^V_(*K0)W+(«)£+ 

£±  =  exp(±^[#3  -  (*K0)3]}  ; 


< 
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the  3\  functions  are  obtained  upon  replacing  #  by  ♦Kj/Kq,  replacing 
•KQ  by  ♦  (for  each  i  =  1,2,  3,4),  and,  for  any  argument  P, 

v  /,«-3 *  105r-6 , 45045--9  .  \ 

=  +  T¥*P  *TTm*  +  )  * 

w*  =  \*1'?gp  ±75otp  ~  screw p  +  ) 


This  result  may  be  used  for  calculation  purposes,  keeping  track  of  the 
relative  magnitudes  of  all  four  arguments  and  the  correct  number  of  terms 
which  therefore  must  be  included.  Retaining  only  the  leading  terms  in  the 
expansion,  we  obtain  the  following  limiting  expression  for  the  transmission 
coefficient: 


T  ~  Y  exp(inQL)  exp^-^ 


which,  for  the  special  case  of  a  straight  trapezoidal  prpfile 


(33) 


-  *Kq 


reduces  to  the  simplified,  useful  expression 


T  ~  y  exp(in0L)  exp(--|)  j«3  -  («K1)3J 


Y  exp(inQL)  exp 


4in0z 


i(,+T^)] 


( 


(34) 


The  convergence  of  (32)  it  such  that  Eqs.  (33)  and  (34)  can  be  used  to  compute 
transmission  when  the  magnitude  of  the  arguments  is  greater  than  two.  The 
corresponding  values  of  L/Kq  would  depend  on  the  conditions  involved  and 
may  be  as  small  as  one  in  a  number  of  cases  of  interest.  Although  (34)  is 
particularly  useful  for  this  purpose  because  of  its  relative  simplicity,  it  is 
limited  to  profiles  which  are  approximately  trapezoidal.  Additional  flexi¬ 
bility  is  provided  by  (33)  which  includes  an  arbitrary  kink  in  the  linear 
representation  of  the  inhomogeneous  region.  In  both  cases,  the  accuracy 
of  the  calculation  can  be  improved,  or  a  smaller  value  of  L/X.^  can  be  con¬ 
sidered  for  a  given  set  of  conditions,  by  including  additional  terms  from  the 
general  expansion  (32). 

A  homogeneous  plasma  slab  geometry  was  used  in  Ref.  2  to  estimate 
the  increase  in  transmission  resulting  from  an  applied  magnetic  field  for 
several  exemplary  re-entry  situations.  One  case  considered  was  that  of  a 
10  degree  cone  at  zero  angle  of  attack  for  a  variety  of  flight  conditions.  In 
the  present  work,  calculations  have  been  made  for  the  same  conditions  using 
the  inhomogeneous  profile  shown  in  Fig.  3.  This  is  a  special  case  of  the 
general  problem  that  is  treated  in  this  section,  which  can  be  shown  to  be  a 
reasonable  approximation  to  the  present  flight  application.  A  comparison 
between  the  more  realistic  inhomogeneous  plasma  values  and  the  previously 
obtained  results  for  a  homogeneous  plasma  is  shown  in  Figs.  4a  through  4c. 

To  discuss  the  effect  of  a  magnetic  field  in  eliminating  "blackout,"  we  must 
establish  a  reference  level  of  acceptable  transmission.  Since  this  value  is 
often  difficult  to  ascertain,  even  for  current  applications  which  are  well 
along  in  the  development  phase,  it  will  be  necessary  to  base  our  remarks  on 
several  possible  assumed  values. 

Consider  first  the  telemetry  frequency  (f  =  240  Me)  and  lower  re-entry 
velocity  given  in  Fig.  4b.  If  10 -db  degradation  in  transmitted  energy  is  ac¬ 
ceptable,  no  magnetic  field  is  required  because  the  possibly  marginal  max¬ 
imum  homogeneous  value  of  9.  75  db  at  150,000  ft  (Bq  =  0)  is  now  replaced 
by  the  corresponding  inhomogeneous  value  of  6.  1  db.  If,  however,  a 
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maximum  of  5  db  it  allowed,  than  blackout  should  occur  at  150,000  ft. 
Moreover,  although  a  magnetic  field  strength  of  250  gauss  would  have  to  be 
applied  to  reduce  the  homogeneous  prediction  to  5  db,  only  50  gauss  is  re* 
quired  to  attain  this  figure  using  the  more  realistic  inhomogeneous  plasma 
calculations.  The  more  severe  re-entry  condition  with  u  =  26,000  fps  is 
shown  in  Fig.  4a (f  *  240  Me).  Although  the  150, 000 -ft  condition  is  still 
more  critical  than  100,000  ft  with  no  magnetic  field,  we  should  note  the 
particularly  strong  dependence  of  db^.  on  the  applied  magnetic  field  in  the 
former  case.  Referring  now  to  inhomogeneous  plasma  results  alone,  we 
see  that  there  is  no  transmission  problem  at  100,000  ft  if  15  db  is  acceptable 
250  gauss  is  needed  to  eliminate  blackout  at  150,000  ft.  For  10  db,  650  gaus 
would  be  required  at  150,000  ft,  while  1350  gauss  would  be  needed  at  100,00( 
ft.  At  the  5  db-level,  the  sharply  reduced  effect  of  Bq  for  100,000  ft  makes 
this  condition  even  more  critical.  Indeed  4200  gauss  would  be  needed  to 
avoid  blackout  in  the  former,  while  1300  gauss  would  be  enough  in  the  latter 
case,  A  qualitatively  similar  situation  is  shown  in  Fig.  4c  when 
f  =  3  X  109  cps  (u  =  26,000  fps). 

The  assumption  of  a  homogeneous  plasma  slab  is  severely  pessimistic 
for  the  re-entry  cone  application.  It  is  of  interest  that  transmission  i3  due 
primarily  to  a  skin  depth  effect.  In  particular,  for  f  =  240  Me  and 
u  =  26,000  fps,  the  fact  that  the  plasma  is  highly  overdense  (Fig.  3)  is  offset 
by  the  correspondingly  small  values  of  L/\q  involved.  The  determination  of 
transmission  for  each  flight  condition  is  greatly  dependent  on  this  effect^ 
whether  the  plasma  is  homogeneous  or  inhomogeneous.  The  influence  of  the 
magnetic  field  on  transmission  is  also  dependent  on  this  situation  since  it  is 
generally  more  pronounced  when  the  level  of  X  is  lower.  *  This  accounts  for 
the  abrupt  variation  of  db^  with  Bq  at  150,000  ft  (Fig.  4a) in  contrast  with  the 
sharply  reduced  effect  at  50,000  ft  where  is  an  order  of  magnitude  larger. 
The  combined  effect  of  all  of  these  quantities  is  shown  in  Figs.  4a  through  4c 
and  the  results  are  discussed  in  the  preceding  paragraph  for  both  the  homo¬ 
geneous  and  more  realistic  inhomogeneous  plasma  models. 
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B .  Exponential-Homogeneous -Exponential  Electron  Density  Distribution 

The  analysis  of  the  preceding  boundary  value  problem  provided  a  clear 
indication  of  the  significant  influence  of  plasma  inhomogeneities  on  propaga¬ 
tion  for  a  range  of  conditions  of  interest.  Detailed  numerical  results  can  be 
obtained  for  specific  problems  also,  approximating  the  actual  electron  density 
distribution  by  the  assumed  profile  (Fig.  1)»  as  was  done  in  the  evaluation  of 
a  re-entering  cone.  It  is  desirable  for  several  reasons  to  consider  yet 
another  formal  inhomogeneous  plasma  slab  problem.  The  need  to  extend  the 
range  of  conditions  previously  considered  to  include  larger  values  of  and 
>as  already  been  noted.  The  present  results  are  also  of  importance  in  the 
small  L/Xq  analysis  of  Section  IV.  Therefore,  let  us  examine  the  formal 
solution  of  Eq.  (11)  for  the  following  assumed  profile  (see  Fig.  5): 


X(z)  =  0 


z  <  0 


X  (1  -  e~az) 
oo' 


a  >  0 


0  <  z  < 


(i-.  *)  .  .,<*< 


x„ 1  -  * 
CD 


(z,  =  L  -  z.) 


=  xoo[l  '  e'a(L'Z>]  .  *2<*<L 

=  0  ,  z  >  L 


(35) 


The  solutions  in  the  several  regions  are  as  follows. 
For  z  <  0: 


F(z) 


„  c<'>  «<2>  .l”0, 


For  0  <  z  <  Zj  ,  Eq.  (1 1)  becomes 
F»  +  nj[l  - 


mXJl  -  e'az)]F  =  0 


(36) 
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Introducing  the  change  of  independent  variables 


2n„ 


1/2  -az/2 


*  =  -rWoo>  '  « 


where  m*^2  is  chosen  so  that  its  imaginary  part  is  positive,  we  obtain 


(37) 


Bessel's  equation 


dZF 

dlZ 


ldF+/ .  a2V.o  a-2'"0 

R'l'Tr  '  ~ 


(1  -  mXJ1^  (38a) 


and,  therefore,  a  general  solution  (a  +  integer)  in  terms  of  Bessel  functions 


F(l)  =  c(3)Ja(I)  +  cWJ  a(f) 


(38b) 


For  Zj  <  z  <  z^: 


az^jl/2 
-az.\T  1/2 


F(z)  =  c<5>  exp(in0)  ft  -  mX^l  -  e  j 

+  c<6>  exp(-in0)£l  -  -  e 

while  for  z^  <  z  <  L  and  z  >  L»,  respectively, 

Fm)  =  c(7)j  mi  +  c(8)j  m 


(39) 


(40) 


t,  =_2(mXoo)1/2  exp[-|(L  -  z)] 


(41) 


F(z)  =  c<9>  e*in°Z 


(42) 
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Continuity  of  F(z)  and  F’(z)  at  z  =  0,  z^,  and  L  gives  rise  to  eight 
equations  for  the  nine  amplitudes  and,  hence,  the  two  ratios  of  interest  (for 
e  ach  wave) 


R 


_  (D"M+  -  C'N+)(C+N+  -  D+M+)  -  y2(D'M"  -  C*N')(C+N'  -  D+M') 

“  1  1  1  * . . . . . . .  '  M . . . . . - . 


(D'M+  -  C"N+)2 


y^d'm' 


C'N')2 


T  = 


-44;Y^a2/-ir2nQ  ^  exp(inQL)  si 


2 

in  a  it 


(D'M+ 


-  Tz - T  -  - 

C  N)4  -  v  (D  M 


C'N')2 


(43) 


(44) 


C*  =  i(mXa>)1/2j;(«)  ±  JQ(f) 

D*  =  l(mXfl0)1/2J^#W  ±  J_Q(*) 

M*  =  (mXJ1^  exp('-Z^)Ja[^  exp( — T-)]  *  exp(‘  ~r)] 

N±  =  (mXoo)1/Z  exp(-“2r)Jltt[^  exp('“z^)]  +  exp(-  ~r~ ) j  • 


2n0  ,  __  .1/2  ,  ..  v  .1/2  2in0,,  v  .  1/2 

0=__(mXQo)  ,  +s(l  -mXj)  ,  Q  =  — — — (1  *  mX(n) 


Y  =  ex] 


p[“o+L(i  -tt)]  ■  xi  =  x„(‘ -  e'*Z') 


In  the  limit  as  a -►  co,  Eqs.  (43)  and  (44)  reduce  to  the  well-known  homogene¬ 
ous  plasma  slab  results. 

The  quantities  and  [see  Eqs.  (29)  and  (30)]  were  computed, 
parametrically,  from  the  solutions  (43)  and  (44).  In  a  similar  manner  to 
the  presentation  of  the  kinked-trapezoid  results,  the  ratio  db,p/(db^,)Q  is 
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plotted  versus  z^/ L  in  Figs.  6a  through  6g  for  a  range  of  conditions.  Four 
values  of  aVg  (characterising  the  profile  shape)  were  assumed,  namely, 
a\g  =  2,4,6,12,  and  the  range  of  and  Y  was  extended  to  include  the  values 
400  and  100,  respectively.  Although  it  was  not  the  purpose  of  this  calcula¬ 
tion  to  compare  the  two  general  distributions  which  have  been  assumed,  such 
a  comparison  was  made  for  several  typical  cases.  As  expected,  the  results 
are  the  same  when  the  profiles  are  ma^  approximately  coincident.  It  is 
also  of  interest  that  the  results  for  the  exponential  geometry  approach  the 
trapezoid  values  when  a  is  sufficiently  small,  as  is  illustrated  graphically  in 
Fig.  6a. 

A  wide  range  of  detailed  profiles  (exponential  and  kinked -trapezoidal 
in  character)  and  values  of  X^,  Z,  and  Y  has  now  been  examined,  parametri¬ 
cally,  for  L/XQ  =  0.  5  and,  to  a  lesser  extent,  L/\q  =  1.  The  resultant 
series  of  illustrative  graphs  which  are  included  in  this  report  can  be  used  to 
effectively  estimate  transmission  for  particular  problems  of  interest,  the 
principal  restriction  in  coverage  being  the  two  values  of  L/Xq  assumed. 
Approximate  calculations  can  be  made  readily,  however,  for  larger  values 
of  WXq  using  the  simplified  Eq.  (34)  or,  if  necessary,  (33)  or  (32).  In  the 
next  section,  a  basic  and  particularly  useful  result  is  established  for  the  small 
L/Xq  regime.  The  present  consideration  of  a  second  class  of  profiles  pro¬ 
vided  the  additional  numerical  foundation  of  an  exact  analysis  which  was 
essential  in  the  verification  of  this  result. 
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SECTION  IV 

PROPAGATION  ACROSS  THIN  INHOMOGENEOUS 
PLASMA  SLABS:*  SMALL  L/XQ 

One  of  the  more  apparent  deficiencies  in  the  extent  of  quantitative  infor¬ 
mation  presented  in  Section  III  lies  in  the  limited  consideration  given  to  small 
values  of  L/\q.  In  the  present  section,  we  will  examine  in  detail  the  propa¬ 
gation  problem  when  L/X^  «  1.  The  desire  to  merely  extend  the  coverage  in 
numerical  information  which  could  be  partially  accomplished  using  a  small  ar¬ 
gument  expansion  of  the  exact  solutions  led  to  the  verification  of  a  rather  basic 
and  particularly  useful  result.  In  this  regard,  it  is  of  interest  to  recall  our 
previously  mentioned  treatment  of  the  re-entry  cone  problem.  All  of  the  con¬ 
ditions  shown  in  Fig.  3  are  in  the  small  L/Xn  range,  even  for  the  higher  fre- 
9  U 

quency  of  f  =  3X10  cps.  Indeed,  it  has  been  observed  that  this  so-called 
skin  depth  effect  is  what  makes  transmission  possible  even  when  the  plasma 
is  highly  over  dense. 

The  present  discussion  will  be  based  primarily  upon  the  use  of  the  exact 
solutions  derived  in  the  preceding  section.  The  following  limiting  situation  is 
offered  to  advance  the  necessary  initial  insight.  It  should  be  emphasized  that 
this  calculation  is  intended  to  serve  as  a  guide  in  the  proper  use  of  our  exact 
solutions  of  the  full  wave  equations;  no  formal  investigation  of  the  problem  is 
intended  along  these  lines. 

Consider  propagation  across  a  vanishingly  thin  plasma  sheath  separating 
two  free  space  regions  as  described  in  Section  II.  In  the  free  space  regions, 
z  <  0  and  z  >  L,  respectively, 


Fj  =0^  °%c<2>ein°2 


(45) 


jQt 

The  author  would  li^e  to  acknowledge  his  appreciation  to  Dr.  Melvin 
Epstein  for  many  fruitful  discussions  concerning  this  subject. 
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(46) 


(3)  -inOz 


F3  =  cw'e 


The  two  waves  propagating  in  the  plasma  (0  <  z  <  L)  satisfy  the  equation 


+  nj[l  -  mX(*)]F2  =  0 


(47) 


We  would  like  to  consider  the  limiting  case  of  L  0,  X  -*oo  (i.e. ,  the 
maximum  value  in  the  layer).  Let  us,  therefore,  represent  the  function 
X(z)  by  a  Dirac  delta  function  at  z  =  0,  6(z), 


X(z)  =  (3fL)6(z) 


IfL  =  /  X(z)dz 

■'o 


(48) 


Substituting  this  into  (47),  we  obtain 

F2  +  n0Fz  •  n0m^L6(z)F2  =  0  .  (49) 


Integrating  from  zero  to  z,  we  see  that 


F2  =  <F2>0+n0m*L<F2>0 


(50) 


2 

where  the  contribution  of  n^  F^  in  (49)  is  zero  in  the  present  limit  since  F^  is 
assumed  to  be  continuous,  and  the  subscript  zero  denotes  the  value  at  z  =  0, 
For  L  -*■<),  therefore,  Fi,  (or  the  magnetic  field)  has  a  simple  discontinuity, 
and  when  the  boundary  conditions  are  applied  we  obtain  the  following  results: 


-  mXn^  L 
2i  +  n^  L 


(51) 
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T  =  1  +  R  = - - - 

1  -  (in0L/2)mX 


«F2»0  _<F2>L  „  'Fi»L  .  _  <Fk>0  .  „ 

— m — ^rr  ■ T  ■  ““oT  •  "^ir  “V1  -  R| 


Evaluating  the  Poynting  vector,  we  obtain 


(52) 


(53) 


*  T 


-k  0Trl2*  |T(|Z) 


M2  +  Z2 

r 


2  2 
uq  +  z 


(54) 


M'  +  [Z  +  (ttLA0)Y]*  M*  +  [Z  +  (irh/\0)Y] 


where  =  1  -  Y,  =  1  +  Y.  (Two  related  problems  of  considerable 
interest  are  briefly  analyzed  in  the  Appendix  using  a  similar  approach.  ) 

We  should  like  now  to  suggest  that  the  problem  of  propagation  across  a 
plasma  slab  in  which  the  ratio  L/\q  is  small,  X(z)  is  arbitrary,  and  the 
maximum  value  of  X  is  large,  is  not  dependent  upon  the  details  of  the  electron 
density  distribution  but  only  on  the  integrated  value  across  the  slab.  Although 
this  assumption  is  made  in  the  analysis  of  radio  waves  from  meteor  trails, 

3 

for  example,  and  in  a  number  of  related  problems  in  classical  physics,  it 
is  our  intention  to  verify  it  using  the  exact  solutions  previously  obtained  in 
Section  III.  As  a  consequence,  a  considerable  simplification  in  the  analysis 
of  specific  propagation  problems  of  this  type  is  obtained. 

An  obvious  deduction  from  the  preceding  remarks  is  that  an  approximate 
solution  to  the  aforementioned  problem  can  be  obtained  from  the  consideration 
of  a  homogeneous  slab  having  the  same  width  L  and  an  equivalent  value  of  TK. 
determined  from  the  relation 
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where  X(z)  ia  the  known  distribution.  For  the  exponential  profile  treeted  in 
Section  III.  B,  therefore, 

^  s  1  '  2^[  T« pUej)  +I^]  •  (56) 

and,  for  the  kinked-trapezoid, 


Therefore,  the  transmission  coefficient  can  be  determined  using  the  corre-  | 

1/2  ' 

sponding  equivalent  index  of  refraction,  n*  =  (1  -  m7)  '  ! 

4K*  exp(innL) 

T  = - « - - - 2 -  (58) 

(IT*  +  1)  exp(inQn*L)  -  (n*  -  1)  exp(-ingn*L) 

It  is  interesting  to  note  that  Eq.  (52)  follows  from  (58)  when  the  magnitude  of 

in^TT^L  is  small  (retaining  only  the  leading  term  in  the  expansion  and  using  the 

assumption  of  small  L/\q,  large  X,  to  further  simplify).  The  basic  premise 

that  propagation  is  principally  dependent  on  the  integrated  value  of  X(z)  when 

L/X.q  «  1  will  now  be  examined  by  comparing  the  results  of  the  approximate 

and  the  exact  solutions.  This  conclusion  is  of  considerable  importance  in 

many  propagation  studies.  The  practical  value  of  using  (52)  or  even  (58)  to  j 

compute  transmission  (or  reflection,  from  the  related  expressions)  for  an 

arbitrary  X(z)  in  the  present  problem  is  evident.  The  following  correlation 

will  also  be  used,  therefore,  to  determine  the  range  of  conditions  for  which  i 

I 

approximate  calculations  can  be  made  on  this  basis.  ] 
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The  transmission  coefficient  has  been  computed  from  the  exact 
solution  (44)  for  the  three  typical  exponential  profiles  shown  in  Fig.  7, 
covering  the  following  range  of  conditions: 


LA,, 

z 

Y 

0.  001 

100,  1000,  10000 

10,  100 

0, 

100 

0.  01 

25,  100,  1000 

1,  10,  100 

0, 

10 

0.  1 

2.  25,  25,  100,  400 

o 

o 

o 

0, 

5,  100 

An  equivalent  homogeneous  plasma  was  determined  using  Eq.  (56)  and  T  was 
computed  from  (58).  The  equivalent  homogeneous  plasma  results  for 
L/\q  =  0.  001  and  0.  01  are  virtually  identical  to  the  exact  values,  indeed 
<dbT)exacl/(dbTbr  =  1. 00,  for  all  cases  considered.  Complete  agreement 
gives  way  to  more  than  acceptable  correlation  for  most  applications  when 
L/Xq  is  increased  to  0.1.  All  three  profiles  and  a  wide  range  of  conditions 
were  included  in  these  calculations  in  order  to  examine  the  applicability  of 
our  approximate  analysis.  Less  than  three-percent  error  was  found  in  the 
computations  involving  profile  I.  Although  the  homogeneous  plasma  values 
‘  of  db,p  for  X  =  Xj  are  already  close  to  the  exact  results  (within  10  percent), 
as  would  be  expected  for  this  profile,  the  additional  accuracy  obtained  using 
X  =  X  is  noteworthy  in  many  applications.  The  discrepancy  increased  in 
general,  to  a  maximum  of  five  percent,  for  profile  II  which  differs  even 
further  from  the  constant  X  =  Xj  distribution.  The  maximum  discrepancy 
noted  in  the  36  cases  considered  for  profile  III  (which  is  now  substantially 
different  than  X  =  Xj)  was  only  seven  percent.  In  general,  even  greater 
accuracy  is  obtained,  for  each  assumed  profile,  as  X^  is  decreased,  for 
fixed  Z,  or  as  Z  is  increased,  for  fixed  Xj. 

In  the  preceding  correlation,  we  established  the  basic  premise  that  the 
propagation  is  dependent  on  If  and  began  to  determine  limiting  conditions 
under  which  the  resultant  convenient  calculation  procedure  can  be  used. 
Consistent  with  the  approach  outlined  at  the  outset  of  this  section,  we  will 


-26- 


not  attempt  to  derive  a  formal  limiting  criterion.  Although  this  would  be 
desirable  and  may  be  examined  using  several  different  expansion  procedures 
suggested  by  the  problem,  the  present  purpose  is  best  accomplished  by  con¬ 
tinuing  to  use  the  same  procedure.  The  exact  solutions,  therefore,  will  now 
be  used  to  provide  the  required  empirical  limiting  information. 

Only  a  relatively  few  calculations  were  made  for  L/Kq  =  0.  5  using  the 
exponential  profile,  in  particular  Xj  =  1. 21,  2.  25,  25;  Z  =  0.  01,  1;  and 

Y  *  0.  It  is  first  apparent  (for  the  extreme  profile  III)  that  agreement  is 
poor,  in  percentage,  when  Z  =  0.  01 ,  Xj  =  1 .  21 .  In  addition,  although  a 
comparison  with  the  corresponding  L/X.q  =  0.  1  results  shows  an  expected 
decrease  in  accuracy,  the  maximum  error  for  the  remaining  conditions  of 
interest  is  still  only  seven  percent.  It  will  be  convenient  to  use  the  analysis 

of  Section  III.  A  to  examine  these  two  points  further.  Transmission  coefficients 
were  obtained  from  the  exact  solution  (26)  for  the  range  of  kinked- trapezoid 
profiles  shown  in  Fig.  8  and  the  conditions 

X1  =  1.44,  2.  25,  4;  Z  =  0.  01,  1;  Y  =  0,  0.5,  1.5,  3 

Using  Eq.  (57),  we  obtain  for  the  equivalent  homogeneous  plasma  slab, 
"X/Xj  =  0.  5(1  +  Xq/Xj  -  Zq/zj).  The  results  are  tabulated  in  Fig.  9.  For 

Y  =  0,  Z  =  1  less  than  three-percent  error  was  obtained  (for  all  profiles)  at 
Xj  =  1. 44  and  2.  25,  and  less  than  seven  percent  at  X^  =  4.  By  contrast,  the 
correlation  experiences  some  difficulties  for  Z  s  0.  01.  The  accuracy  is 
generally  improved  as  Xj  is  increased  (depending  on  the  profile),  particularly 
for  Z  =  0.  01;  indeed  the  error  obtained  at  Xj  =  4,  Z  =  0.  01,  now  ranges  from 
one  to  13  percent  for  all  cases  (except  one).  At  the  same  time,  note  also  that 
for  a  fixed  value  of  X ^  there  is  a  marked  improvement  in  the  correlation  when 
either  Xq/Xj  is  increased  with  Zq/zj  held  fixed  or  Zq/zj  is  decreased  with 
Xq/Xj  fixed.  The  principal  factor  involved  appears  to  be  the  abrupt  depend¬ 
ence  of  €  j  on  X  near  the  resonant  condition  when  Z  is  small.  Part  of  the 
plasma  is  always  overdense  because  we  have  assumed  Xj  >  1 .  A  smaller 
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value  of  Xq/Xj  and/or  a  larger  value  of  will  reduce  the  ratio  7/Xj, 

the  minimum  value  being  0.  2  for  the  profiles  considered.  In  this  case,  for 
Xj  =  1.44,  7=  0.  288  and  the  equivalent  homogeneous  plasma  will  not  include 
any  resonant  effect.  Although  the  percentage  error  is  large  in  this  case,  it 
is  interesting  to  observe  that  the  magnitude  of  db,p  is  quite  small  which 
therefore  may  be  more  of  academic  than  practical  concern. 

At  the  other  extreme,  7=1.  152,  the  resonant  condition  is  more 
accurately  taken  into  account,  and  an  error  of  only  12  percent  is  obtained. 

Note  the  magnitude  of  db^,  involved.  As  X j  is  increased,  the  region  in  which 
X  is  of  order  one  represents  a  smaller  portion  of  the  total  distribution.  Its 
effect  on  transmission  will  be  reduced,  and  a  qualitatively  similar  but 
substantially  improved  correlation  is  shown  for  =  2.  25.  As  was  previously 
noted,  the  over-all  agreement  is  quite  good  when  X^  =  4.  The  major  effect 
described  above  is  still  evident  since  for  the  extreme  profile  when  7  =  0.  8  the 
percentage  error  is  sizable,  whereas  good  agreement  is  obtained  in  all  of  the 
remaining  cases  since  the  equivalent  value  of  7  is  further  from  the  resonant 
value.  The  maximum  discrepancy  is  less  than  seven  percent  at  Z  =  1  for  all 
cases  including  a  magnetic  field.  The  way  in  which  X^  and  Y  combine  to 
effect  the  correlation  may  be  seen  from  the  detailed  results  given  in  Fig.  9. 

The  presence  of  a  magnetic  field  modifies  the  preceding  discussion  for 
2  =  0.01  in  certain  instances.  Although  the  qualitative  conclusions  often  apply, 
it  is  best  to  refer  to  the  detailed  correlation  to  evaluate  the  prospective 
accuracy  of  the  approximate  procedure  for  a  specific  combination  of  Xj  and  Y 
of  interest. 

The  following  conclusions  can  be  made  in  summary.  Using  the  exact 
solutions  derived  earlier  in  this  report,  we  have  shown  that  the  transmission 
across  a  thin  (L/\q  «  1)  inhomogeneous  plasma  sheath  depends  on  the  total 
integrated  value  of  electron  density,  not  on  the  detailed  distribution.  Such  an 
assumption  is  often  the  essential  starting  point  for  investigations  of  propagation 
in  this  limit  and  is  of  particular  value  in  the  present  analysis.  Accurate 
transmission  calculations  can  be  made  on  this  basis  for  L/Kq  <  0.  1  using  the 
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considerably  simplified  homogeneous  plasma  slab  expression  (58)  and  an 
equivalent  constant  value  of  X  determined  from  (55).  Indeed  for  the  broad 
range  of  conditions  considered  in  establishing  this  fact,  there  is  virtually 
complete  agreement  between  the  exact  and  the  approximate  results  for 
WKq  *  0.001  and  0.01.  When  L/XQ  was  increased  to  0.  1,  the  maximum 
discrepancy  found  in  the  computations  for  108  different  cases  (including 
three  profiles),  covering  the  full  spectrum  of  transmission  levels  of  possible 
interest,  was  only  seven  percent.  In  general,  even  greater  accuracy  is 
obtained  when  the  profile  approaches  the  constant  X  =  X^  distribution  or,  for 
an  assumed  profile,  when  either  X^  is  decreased  or  Z  is  increased.  For 
L/Kq  =  0.  5,  smaller  values  of  X^,  Z,  and  Y  were  used,  along  with  a  number 
of  profiles,  and  a  maximum  error  of  seven  percent  was  obtained  when  Z  =  1. 
The  correlation  at  Z  =  0.  01  was  substantially  affected  by  the  values  of  Xj 
and  Y  and  by  the  profile.  Although  the  previous  discussion  for  Y  =  0  is 
informative,  it  is  best  to  refer  to  the  detailed  results  of  Fig.  9  to  evaluate 
a  particular  condition  of  interest. 

It  should  be  remarked  that  we  have  not  attempted  rigorously  to  determine 
a  formal  limiting  condition  for  the  approximate  calculation  procedure.  Although 
its  accuracy  (better  than  93  percent)  was  established  in  detail,  in  an  empirical 
manner  for  virtually  any  condition  of  practical  interest,  even  when  L/X^  =  0.  1, 
only  limited  calculations  were  made  for  larger  values  of  L/  X^  where  the 
procedure  will  begin  to  fail.  However,  for  L/Xq  =  0.  5  (and,  to  a  lesser  extent, 
Lj/\q  =  1)  transmission  can  be  adequately  determined  from  the  detailed  con¬ 
siderations  of  Section  III  which  includes  a  large  measure  of  arbitrariness  in 
the  assumed  distribution  of  X(z)  and  *  broad  parametric  treatment  of  the 
remaining  quantities.  The  evaluation  of  our  appproximate  procedure  at 
L/Xq  =  0.  5  is  intended  to  complement  these  results  when  calculations  are 
required  for  intermediate  values  of  the  parameters  or  for  transitional  values 
of  L/Xq,  say  between  0.  1  and  1.0.  We  shall  now  illustrate  the  way  in  which 
the  two  procedures  can  be  used  to  complement  one  another.  Recall  that  accu¬ 
rate  results  (less  than  seven- percent  error)  were  obtained  from  the  approxi¬ 
mate  calculations  at  L»/Xq  =  0.5  when  Z  =  1.  Thus,  transmission  results  can  be 


-29- 


( 


computed  at  Z  =  1  for  the  intermediate  valuee  of  which  are  not  included 
in  the  graphs  obtained  in  Section  III.  A.  Estimates  also  can  be  made  on  this 
basis  for  values  of  Z  which  are  close  to  one,  covering  the  entire  range  of 
considered.  On  the  other  hand,  the  dependence  on  Z  shown  in  these  graphs 
can  be  used  to  estimate  transmission  when  Z  is  small,  aided  to  some  extent 
by  an  approximate  calculation  whose  accuracy  can  be  evaluated  from  Fig.  9. 

A  vivid  demonstration  of  the  obvious  utility  of  the  simplified  small 
L/Xq  expressions  follows.  A  triangular  distribution  of  X(z)  was  previously 
assumed  to  represent  the  plasma  about  a  re-entry  cone  (Fig.  3).  With  the 
aid  of  a  high-speed  computer,  the  exact  solution  derived  in  Section  III.  A  was 
used,  therefore,  to  calculate  transmission.  The  results  obtained  when 
f  =  240  Me  can  be  easily  reproduced  (within  several  percent)  using  the 
limiting  algebraic  expression  (54).  Satisfactory  agreement  was  not  obtained 
at  f  =  3,  000  Me  when  (54)  was  no  longer  appropriate,  i.  e.  ,  for  conditions 
where  the  magnitude  of  in^lT^L  is  not  small.  Excellent  agreement  with  the 
exact  calculations  was  obtained  in  all  cases  using  Eq.  (58).  It  follows  from 
the  general  correlation  of  (58)  with  the  exact  solutions  of  Section  III  that  (54) 
may  be  used  to  simplify  computations  even  further,  when  applicable.  It  is 
perhaps  redundant  to  add  that  homogeneous  plasma  slab  calculations  can  be 
obtained,  in  general,  from  (54)  rather  than  from  the  more  complicated  (58) 
when  |inQn*L|  is  small.  The  result  may  still  be  substantially  less  than 
perfect  transmission  in  this  situation  since  a  small  value  of  may  be 
offset  by  a  correspondingly  large  value  of  |n*|  with  the  product  still  remaining 
less  than  one. 
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SECTION  V 


CONCLUSIONS 

Exact  solutions  for  the  reflection  and  transmission  coefficients  have 
been  obtained  for  several  exemplary  transition  zones  between  dielectric 
and  dissipative  gases.  Normal  incidence  into  a  stratified  plasma  slab  is 
assumed  so  that  the  electromagnetic  waves  are  propagating  parallel  to  the 
free  electron  density  gradients.  The  magnetic  field  is  constant  and  is  applied 
in  the  propagation  direction.  Two  general  electron  density  distributions  are 
considered  in  which  both  the  shape  and  the  extent  of  the  inhomogeneous  profile 
may  be  varied. 

(i)  Kinked- trapezoid  (Fig.  1), 

(ii)  Exponential- homogeneous -exponential  (Fig.  5). 

The  solutions  derived  on  this  basis  are  used  to  evaluate  the  effect  of  an 
inhomogeneous  plasma  on  electromagnetic  wave  propagation,  parametrically, 
both  with  and  without  an  applied  magnetic  field.  At  the  same  time,  previously 
obtained  homogeneous  plasma  calculations  for  a  re-entry  cone  are  refined 
using  the  analysis  of  problem  (i). 

The  kinked-trapezoid  geometry  is  considered  in  Section  III.  A,  and  the 
resultant  expressions  for  the  reflection  and  transmission  coefficients  are 
given  by  Eqs.  (25)  and  (26),  respectively.  Parametric  calculations  were 
made  for  L/Kq  =  0.  5,  1.0  including  moderately  overdense  plasmas,  i.  e.  , 

1  <Xj  £  4,  comparable  values  of  0  <  Y  ^  3  required  to  improve  transmission 
to  the  level  of  practical  interest,  and  0.  01  <  Z  <  1.  Illustrative  results  for 
db^,,  normalized  by  the  corresponding  homogeneous  value  are  shown  in 
Figs.  2a  through  2k  varying  the  details  of  the  assumed  profile. 

The  second  general  profile  was  examined  in  Section  III.  B  with  the 
solutions  given  by  Eqs.  (43)  and  (44).  The  resultant  calculations  for 
L/Kq  *  0.  5,  1.0,  including  larger  values  of  Xj  and  Y,  are  illustrated  in  a 
similar  manner  in  Figs.  6a  through  6g.  The  wide  range  of  detailed  profiles 
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(exponential  and  kinked- trapezoidal  in  character)  and  values  ofXj,  Z,  Y 
covered  in  the  two  sets  of  graphs  can  be  used  to  estimate  transmission,  the 
principal  restriction  being  the  two  values  of  L/Xq  assumed. 

Transmission  calculations  for  larger  values  of  L/X^  can  be  directly 
obtained  from  the  asymptotic  expansion  of  Eq.  (26).  The  resultant  limiting 
expressions  (34)  and  (33)  are  applicable  for  values  of  L/X^  as  small  as  one, 
depending  on  the  conditions  involved,  due  to  the  rapid  convergence  of  (32)  when 
the  magnitude  of  the  arguments  is  of  the  order  of  two.  Although  (34)  is 
particularly  useful  for  this  purpose  because  of  its  relative  simplicity,  it  is 
limited  to  profiles  which  are  approximately  trapezoidal;  additional  flexibility 
is  provided  by  (33)  which  includes  a  kink  in  the  profile.  In  both  cases  the 
accuracy  of  the  calculation  can  be  improved,  or  a  smaller  value  of  L/Xq  can 
be  considered  for  a  given  set  of  conditions,  by  including  additional  terms 
from  the  general  expansion  (32). 

The  detailed  exact  solutions  of  problems  (i)  and  (ii)  were  used  for  yet 
another  purpose  in  the  analysis  of  Section  IV.  In  that  section  the  problem 
of  propagation  across  thin  inhomogeneous  plasma  slabs  \\  is  treated. 

Numerical  calculations  using  these  exact  solutions  were  the  principal  means 
of  verifying  the  resultant  small  L/Xq  analysis.  A  brief  consideration  of  the 
limiting  problem  of  a  vanishingly  thin  plasma  sheath  (L  -*0,  X  -*ao)  resulted 
in  the  solutions  (51)  and  (52)  for  the  reflection  and  transmission  coefficients 
and  (54)  for  the  transmitted  energy.  (Two  related  problems  of  interest  are 
considered  in  the  Appendix.  )  On  this  basis,  it  was  assumed  that  when  L/Xq 
is  small  and  X(z)  is  arbitrary,  propagation  is  not  dependent  upon  the  details 
of  the  electron  density  distribution  but  only  on  the  integrated  value  across  the 
slab.  A  considerable  simplification  in  the  analysis  of  problems  of  this  type 
should  be  obtained  therefore  from  the  solution  of  an  equivalent  homogeneous 
slab  problem,  Eq.  (58),  with  X  =  X  =  constant  defined  in  (55).  The  validity 
of  these  assertions  was  examined  by  comparing  the  appropriate  exact  and 
approximate  expressions  for  a  wide  range  of  conditions.  The  virtually 
identical  results  obtained  for  L/X^  «  1  (in  particular,  L»/Xq  =  0.001  and  0.01) 


-32- 


established  the  basic  premise  regarding  the  dependence  on  X  which,  it  may 
be  noted,  can  be  used  in  a  wide  class  of  problems  apart  from  the  present 
general  boundary  value  problem  outlined  in  Section  II, 

Complete  agreement  became  more  than  acceptable  correlation  when 
L/\Q  was  increased  to  0.  1,  Indeed,  a  maximum  discrepancy  of  seven 
percent  was  found  in  the  calculations  covering  the  full  range  of  transmission 
levels  of  possible  interest  corresponding  to  the  conditions,  =  2.  25,  25, 

100,  400;  Z  =  0.  01,  1,  10;  Y  =  0,  5,  100.  A  further  correlation  at 
L/\Q  =  0.  5  was  derived  in  which  the  degree  of  accuracy  depended  on  the 
values  of  these  parameters.  This  information  is  of  value  as  a  complement 
to  the  extensive  results  given  in  Section  III  for  L/\q  =  0.  5  (and,  to  a  lesser 
extent,  L/X.^  =  1)  when  calculations  are  required  at  other  values  of  X^,  Z,  Y, 
or  in  the  transitional  region  of  0.  1  to  1 . 0. 

2 

The  problem  of  a  re-entry  cone  which  was  previously  treated  assuming 
a  homogeneous  plasma  is  considered  once  again  in  Section  III.  A,  using  a 
triangular  distribution  of  X(z).  The  geometry  and  cases  are  specified  in 
Fig.  3  and  the  results  are  shown  in  Figs. 4a  through  4c.  The  increase  in 
transmission  predicted  by  the  more  realistic  inhomogeneous  plasma  formu¬ 
lation  is  quite  pronounced.  Further,  the  magnetic  field  required  to  provide 
a  specified  acceptable  level  of  transmission  can  be  substantially  larger 
assuming  a  homogeneous  plasma.  In  Fig.  4af  for  example  (f  =  240  Me, 
u  =  26,  000  fps),  at  100,  000  ft,  1200  gauss,  3600  gauss,  and  7800  gauss  would 
be  required  to  obtain  15  db,  10  db,  and  5  db,  respectively,  for  a  homogeneous 
plasma.  The  corresponding  values  specified  by  the  inhomogeneous  plasma 
calculations  are  0,  1400,  and  4200  gauss,  a  significant  decrease  of  practical 
interest.  The  detailed  values  for  this  typical  re-entry  application  are  shown 
in  the  aforementioned  graphs.  The  obvious  utility  of  the  results  of  the  small 
L/\q  analysis  is  clearly  demonstrated  by  its  application  to  the  re-entry  cone 
problem.  Good  agreement  was  obtained  between  the  exact  values  shown  in 
Figs. 4a  and  4b  and  the  corresponding  limiting  algebraic  expression  results 
[Eq.  (54)]  .  Excellent  agreement  was  obtained  in  all  cases  when  the  more 
general  approximate  expression,  (58),  is  used,  including  the  conditions  in 
Fig.  4c  at  which  |inQn*L|  is  no  longer  small. 
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Fig.  2b.  Transmitted  energy,  kinked  trapezoid  X(z)  (continued). 
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Fig.  2d.  Transmitted  energy,  kinked  trapezoid  X(z)  (continued). 
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Fig.  2e.  Transmitted  er.ergy,  kinked  trapezoid  X(z)  (continued). 
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Fig.  2i.  Transmitted  energy,  kinked  trapezoid  X(z)  (continued). 
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Fig.  2k.  Transmitted  energy,  kinked  trapezoid  X(z)  (continued). 
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Fig.  3.  Re-entry  10°  cone  at  zero  angle  of  attack. 
Inhomogeneous  plasma  profile,  flight  conditions. 


Altitude 
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Fig.  4a.  Transmitted  energy,  re-entry  10°  cone  at  zero  angle  of  attack. 
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Fig.  4b.  Transmitted  energy,  re-entry  10°  cone  at  zero  angle  of  attack. 

(continued) 
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Fig.  6d.  Transmitted  energy;  exponential  X(z)  (continued). 
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I  Fig*  6f.  Transmitted  energy;  exponential  X(z)  (continued). 
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I  Fig.  6g.  Transmitted  energy;  exponential  X(s)  (continued). 
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rig.  8.  Kinked -trapezoid  profiles  used  for  correlation  of  equivalent 
homogeneous  plasma  transmission  calculations  with  exact  values. 
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Fig.  9a.  Correlation  of  equivalent  homogeneous  plasma  transmission 
calculations  with  exact  values  for  a  kinked-trapesoid  profile. 
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Fig.  9b.  Correlation  of  equivalent  homogeneous  plasma  transmission 
calculations  with  exact  values  fora  kinked  -traps  sold  profile  (continued), 
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Fig.  9c.  Correlation  of  equivalent  homogeneous  plasma  transmission 
calculations  with  exact  values  fora  kinked  -tr  ape  acid  profile  (continued). 
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Fig.  9d.  Correlation  of  equivalent  homogeneous  plasma  transmission 
calculations  with  exact  values  for  a  kinked -trapezoid  profile  (continued). 
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APPENDIX 


We  shall  consider  first  the  propagation  of  a  horizontally  polarized  wave 
which  is  obliquely  incident  on  a  vanishingly  thin  isotropic  plasma  sheath.  The 
governing  equation  for  the  electric  field  intensity  is  now  given  by^ 

Ey  +no(C'2  *  Tj)  Ey  =  0  '  C  =  co80I  ,  (A-l) 


where  0j  is  the  angle  of  incidence,  i.  e.  ,  the  angle  formed  by  the  wave  normal 
of  the  incident  wave  and  the  positive  z-axis.  A  similar  procedure  to  the  one 
employed  at  the  outset  of  Section  IV  leads  quite  readily  to  the  following  results: 


-Xn0L/U 

2iC  +  5nQL/U 


(A- 2) 


1  -  (inQL/2C)Y/U 

_ _ 1  +  Z2 _ 

l  +  [z  +  (WC)(LA0)Y]  2 


(A- 3) 


(A- 4) 


These  expressions  reduce  to  the  corresponding  ones  for  the  case  of  normal 
incidence  (0^  =  0)  derived  in  Section  IV. 

The  function  c  a  maximum  at  0j  =  0,  decreasing  monotonicaliy 

to  zero  at  0  j  =  *90  degrees.  Since  c  ^  depends  only  on  X,  we  are  led  again  to 
the  assumption  of  an  equivalent  homogeneous  plasma  slab,  in  which  case, 

(4/C)(K*2  -  S2)1  ^2  exp(inQCL) 

[1  +  (**Z  -  S2)1/i!/C]*  exp[in0MTt*i  -  sV/Z]  -  [1  -  <W*Z  -  SZ)l/2/C]  exp  [-in0L(H*Z  -  S2)l/Z] 


(A-5) 
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where  S  =  sin  6j  and  n*^  =  (1  -  X/U).  This  expression  will  reduce  to  (A-3) 
in  the  present  problem  when  the  magnitude  of  in^rT^L  is  small.  The  extensive 
correlation  made  in  Section  IV  for  0^=0  between  the  values  obtained  from 
the  approximate  expression  (58)  [or  (52),  if  applicable]  and  the  appropriate 
exact  results  can  now  be  «.sed  to  justify  the  use  of  (A- 5)  [or  (A-3))  in  the 
following  way:  For  given  Y,  Z,  and  L/\Q  (Y  =  0),  Eq.  (A-3)  is  the  same  as 
(52)  if  we  consider  a  new  value  of  LAq  which  is  modi^iecl  by  the  factor  1/C. 

It  might  be  anticipated  that  the  excellent  agreement  obtained  in  this  correlation 
for  L/Kq  ^  0.  1  and  the  subsequent  discussion  for  L/X.^  =  0.  5  are  an  indication 
of  the  accuracy  afforded  by  (A- 5)  in  the  present  problem  for  the  corresponding 
conditions 


v —  <  0.  1  cos  0 T  *  t-—  =0.5  cos  0T  .  (A-6) 

*0  1  \)  A 

The  preceding  extension  to  the  case  of  oblique  incidence  (Y  =  0)  is 
certainly  tentative  in  many  respects.  It  is  felt,  however,  that  (A-3)  or  (A- 5) 
can  be  used  to  give  approximate  results  in  many  applications,  with  qualitative 
limitations  provided  by  the  correlation  of  Section  IV  along  with  the  modifica¬ 
tion  (A-6).  As  a  final  footnote,  we  shall  consider  the  following  even  more 
tentative,  but  possibly  useful,  reasoning  when  Z  =  Z(z).  In  general 

w"  +  n2(C2  -  mX>w  =  0  ,  (A-7) 

where  w  =  F  and  C  =  1  (normal  incidence)  when  there  is  an  applied  magnetic 
field;  w  =  for  Y  =  0;  and  m  =  (U  ?  Y)  *  is  now  a  function  of  z.  The 
previous  limiting  calculations  can  be  altered  accordingly  and  we  obtain 

_Ino 

R  =  "ZTCTTn^  ’  <A-8> 
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! 


whe  re 


(A- 9) 


(A-  10) 


It  seems  reasonable  to  suggest  that  these  limiting  results  may  be  used  to  make 

approximate  calculations  including  a  variable  Z(z)  when  L/\q  «  1  and 

|inAL(mX)  I  is  less  than  one. 

'  0  v  'max  1 
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